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Abstract 

- - - , We give a linear-time algorithm that checks for isomorphism between two — 1 matrices 

^ ' that obey the circular-ones property. This algorithm leads to linear-time isomorphism algo- 

f— ^ , rithms for related graph classes, including Helly circular-arc graphs, F circular-arc graphs, 

^SJ ' proper circular-arc graphs and convex-round graphs. 

1 Introduction 
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Two graphs G and G' are isomorphic if there is a bijection vr from the vertex set of G to the 

vertex set of G' that satisfies the following condition: for all n, u S V{G)^ u is adjacent to v 

in G if and only if 'k{u) is adjacent to ^{v) in G' . The graph isomorphism problem consists of 

determining whether two input graphs are isomorphic. The graph isomorphism problem is in 

NP, but it is neither known to be NP-complete nor known to be in P. However, the restriction of 

Q ■ the problem to members of graph classes with certain special topological properties is known to 

result in polynomial or even linear-time algorithms. The subject of this paper is the development 

of such algorithms for a variety of such classes. 

K^ ■ We show linear-time isomorphism algorithms for Helly circular-arc graphs, F circular-arc 

^vj ■ graphs, proper circular-arc graphs and convex-round graphs. The common building block for 

00 ! all of these algorithms, is a linear-time isomorphism algorithm for binary matrices that obey 

'^ ' the circular-ones property, which we show. 

cn . In order to explain our results, we first establish some basic terminology. A matrix is a 

binary matrix if all of its entries are or 1. The adjacency matrix of a simple graph G is a 
binary square matrix M that has 1 in row i, column j, if vertex i is adjacent to vertex j, and 
otherwise. The augmented adjacency matrix of G is the adjacent matrix of the graph with I's 
on the main diagonal. 

A clique matrix of a graph is a binary matrix that has one row for each vertex and one 
cH • column for each maximal clique, and a 1 in row i, column j if vertex i is a member of maximal 

clique j. 

A consecutive- ones matrix is a binary matrix whose columns can be ordered such that, 
in every row, the I's are consecutive. Such an ordering of the columns is a consecutive- ones 
ordering. A circular-ones matrix is a binary matrix whose columns can be ordered such that, 
in every row, either the O's are consecutive or I's are consecutive; equivalently, the I's are 
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consecutive modulo the number of columns, and the block of I's is allowed to "wrap around" 
from the rightmost to the leftmost column. Such an ordering of the columns is a circular-ones 
ordering. It is easily seen that the class of consecutive-ones matrices is a proper subclass of the 
class of circular-ones matrices |13j . 

If a circular-ones ordering of a matrix with n rows is known, then it can be represented in 
0{n) space by recording the number of columns and listing, for each row, the columns of the 
first and last 1 in the row. The column of the first 1 is the column where a 1 follows a zero 
and the column of the last 1 is the one where a 1 is followed by a 0. (The column of the last 1 
precedes that of the first if the row wraps around the end of the matrix.) A row that has only 
O's or only I's can be represented with a suitable code, such as a single or a single 1. Let 
us call this a succinct representation of a circular-ones matrix. There may be many succinct 
representations, since there may be many circular-ones orderings of the matrix. 

The intersection graph of a family of sets has one vertex for each set in the family and an 
edge between two vertices if the corresponding sets intersect. 

A circular-arc graph is the intersection graph of arcs on a circle. If we restrict circular-arc 
graphs to intersection graph of arcs on a circle such that no arc contains another, we get proper 
circular- arc graphs |13j . 

Another subclass of circular- arc graphs is the Helly circular- arc graphs, sometimes called 
6 circular- arc graphs. A circular- arc model has the Helly property if every family of pairwise 
intersecting arcs has a common intersection point. Such a model is called a Helly circular-arc 
model. A circular-arc graph is a Helly circular-arc graph if there exists a Helly circular-arc 
model of it. In [11], it is shown that a graph is a Helly circular-arc graph if and only if the 
clique matrix satisfies the circular-ones property. 

The graphs whose augmented adjacency matrices satisfy the circular-ones property are T 
circular-arc graphs, also called concave-round graphs. This graph class is also a subclass of 
circular-arc graphs [3T]. 

The complement G of an undirected graph G has the same vertex set, and an edge between 
two vertices if and only if there is no edge between them in G. 

The complement of a F circular-arc graph has the circular-ones property for its adjacency 
matrix. This class of graphs is called convex-round graphs \2\ . 

Wu |32j presented the first polynomial algorithm for circular-arc graph isomorphism, but 
later Eschen [9] claimed to find a flaw in it. Hsu [15] presented an 0{nm) isomorphism algorithm 
for circular-arc graphs where n denotes the number of vertices and m denotes the number of 
edges in a graph. In Section [7] we give a counterexample to the correctness of this algorithm. 
We also describe there a suggestion given by Hsu for a possible fix for the algorithm. Therefore, 
there are currently no known efficient isomorphism algorithms for circular-arc graphs. Some 
subclasses of circular- arc graphs do have efficient isomorphism algorithms. Interval graphs [23], 
co-bipartite circular arc graphs [9], and proper circular-arc graphs |22] all have linear-time 
isomorphism algorithms, while F circular-arc graphs have an 0{n'^) isomorphism algorithm [S]. 

Two binary matrices Mi and M2 are isomorphic if there exists a permutation r of the rows 
of Ml and a permutation vr of its columns that makes Mi identical to M2. If vr is known, r 
is trivial to find by matching up identical rows of Mi and M2. Therefore, abusing notation 
somewhat, we will sometimes call vr an isomorphism from Mi to M2, omitting mention of r, 
and treat a matrix as a multiset of row vectors. 

The results of the paper are organized as follows. In Section [21 we give basic definitions 
and review the PC tree of a circular-ones matrix |17t 128] . This gives a representation of all 
circular-ones orderings of a circular-ones matrix. 

In Section [3l we present a notion of quotient labels on the PC tree, which were developed 
in [8]. The matrix can be reconstructed from the tree and its labels, establishing the quotient- 



labeled PC tree as a unique decomposition of a circular-ones matrix. 

In SectionHl we give an algorithm that uses the quotient-labeled PC tree to test isomorphism 
of circular-ones matrices that was also developed in [8j. We define a notion of isomorphism of 
quotient-labeled PC trees, show that two circular-ones matrices are isomorphic if and only if 
their quotient-labeled PC trees are isomorphic, and reduce the problem to testing whether the 
two matrices' PC trees are isomorphic. The running time is linear in the number of rows, 
columns and I's of the matrix if the circular-ones orderings are not provided, or linear in the 
number of rows if a succinct representation of circular-ones matrices is provided. 

In Section [5l we reduce the problem of testing whether two Helly circular-arc graphs are 
isomorphic to testing whether two circular-ones matrices are isomorphic. This gives a bound 
that is proportional to the number of vertices and edges, or just proportional to the number of 
vertices, depending on whether the graphs are represented with adjacency lists or with suitable 
sets of circular arcs. A preliminary version of part of the results of this section appeared in [21] . 

In Sectional we use the fact that testing isomorphism of T circular- arc graphs and of convex- 
round graphs reduces to testing isomorphism of circular-ones matrices, giving an 0{n + m) or 
an 0{n) bound, depending on whether a succinct representation is given. This leads to a new 
algorithm for testing the isomorphism of proper circular arc graphs, which can run in 0{n) time 
if two circular arc models are given. 

In Section [7] we discuss the circular-arc isomorphism algorithm of [T5] . We show a coun- 
terexample for this algorithm, and give a direction suggested by Hsu for a possible fix. 

2 Preliminaries 

We consider simple undirected graphs G and G'. We denote the number of vertices of G by n 
and the number of edges by m. We assume that G' has the same number of vertices and edges 
as G, since otherwise it is trivial to see that the graphs are not isomorphic. 

In this paper, we also let n denote the number of rows of a binary matrix. (Many of the 
matrices we deal with are derived from graphs, and have one row for each vertex of the graph.) 
The size of a binary matrix M, denoted size{M), is the number of rows plus the number of 
columns plus the number of I's; this is proportional to the number of words required to store 
the matrix using a standard sparse-matrix representation. We will say that an algorithm whose 
inputs are binary matrices runs in linear time only if it runs in time that is linear in this measure 
of the size of the matrices. 

By N{v), we denote the set of neighbors of a vertex v. By N[v], we denote the closed 
neighborhood {v} U N{v) of v. If [/ is a set of vertices then N[U] is the union of N[v] for all 

veU. 

An (unrooted) tree T is an undirected graph that is connected and has no cycles. Rooting 
a tree T at node w consists of orienting all edges so that they are directed from vertices that 
are farther from w to vertices that are closer to w, yielding a directed graph. This gives each 
vertex u ^ w a unique outgoing edge {u, v). The neighbors of u can be classified as the (unique) 
parent of u (the unique neighbor that is closer to w) and the children of u. Once a tree has been 
rooted, we will continue to refer to the parent and children of u as its neighbors. 

A circular-arc model of G is a mapping from the vertices of G to arcs on a circle such that two 
vertices are adjacent if and only if the corresponding arcs intersect. We represent a circular-arc 
model ^ by a cyclic doubly-linked list of the endpoints arcs. Each vertex of G has two endpoints 
in the model, one of them is a clockwise endpoint and the other is a counterclockwise endpoint. 

A proper circular-arc model is a circular-arc model in which no arc contains another. An 
interval model is a circular-arc model whose circle has some point that is not contained in any 
arc. 



Let t = {ti,t2, ■ ■ ■ ,tk) he a list, where each tj is a tuple of integers. By sorting a list of tuples, 
we mean that we rearrange the order of the ij's so that they are in increasing lexicographic order. 
Let Li and L2 be two lists of tuples. Let i be the first position in which they differ. Then we say 
that Li precedes L2 lexicographically if the tuple at position i of Li lexicographically precedes 
the tuple at position i of L2, or else Li has length i — 1. 

2.1 Bipartitive decomposition trees 

Let a bipartitive tree T on universe V be an undirected tree such that its leaves are the elements 
of V, all internal nodes have degree at least three, and each internal node is labeled either prime 
or degenerate. Let a neighbor set of a node u denote the set of leaves in a tree of the forest that 
results when u is removed from T. Since all internal nodes of T have degree at least three, all 
members of a neighbor set are leaves of T, so each neighbor set is a subset of V. The set family 
J-{T) represented by T consists of the following sets: 

• A neighbor set of a prime node or the union of all but one of the neighbor sets of a prime 
node; 

• Any union of at least one and fewer than all neighbor sets of a degenerate node. 

Not all set families can be represented by a bipartitive tree. Next, we characterize those 
that can. 

The symmetric difference X/\Y of two sets X and Y is the set {X \ y) U (y \ X). Let 
us say that two subsets X and Y of universe V strongly overlap ifXny, X\Y , Y \X and 
XTJY = V\{X{JY) are all nonempty. 

A bipartitive set family on universe y is a set family F with the following properties: 

• 0,y 07" 

• {x} € J- for all X € y 

• X G J" for all X G J" 

• Whenever X, y G -F strongly overlap, X r\Y , XiJY and Xl^Y are all members of J-. 

Theorem 2.1 ([7J). If T is a bipartitive tree on universe V , then J-(T) is a bipartitive set 
family on universe V. Conversely, if T is a bipartitive set family, then there exists a unique 
bipartitive tree T such that T = J-{T) . 

2.2 PC trees 

Let M be a circular-ones matrix. A row of M can be thought of as the bit- vector representation 
of a set, that is, it is the set X of columns of M where the row has a 1. Let V denote the 
columns of M and let TZ denote the family of sets represented by the rows. Note that 7^ is a set 
family on universe V . Let Af{TZ) denote the family of subsets of V, excluding and V itself, 
that do not strongly overlap any member of 7^. 

Lemma 2.2 (|17]). J\f{TZ) is a bipartitive set family on universe V. 

It follows that M{TZ) is represented by a bipartitive tree. For historical reasons, the prime 
nodes are known as C nodes, the degenerate nodes are known as P nodes and the bipartitive 
tree is called a PC tree \\.7\ 128]. Figure [1] gives an example. When drawing a PC tree, we use 
the convention of representing a C node with a double circle and a P node with a dot. In this 



figure, the neighbor sets of c are {1,2,3}, {4}, {5}, {6}, and {7,8,9,10}. Each of these sets 
and the union of all but any one of these sets is a member of Af{TZ). For example, the neighbor 
set {7, 8, 9, 10} is a subset of rows {3, 5, 13, 14}, it contains rows {6, 7, 8}, it does not strongly 
overlap row 9 because the union of row 9 and {7, 8, 9, 10} is the entire column set, and it is 
disjoint from all other rows. Since it does not strongly overlap any row, it is a member oi MiJZ). 
Similarly, the union of all neighbor sets other than {7, 8, 9, 10} (the complement of {7, 8, 9, 10}), 
is a member of M(JZ). However, because c is a C node (prime) the union of neighbor set {6} 
and {7, 8, 9, 10}, which is neither the union of one neighbor set nor the union of all but one, is 
not a member of M(7Z). This is verified by observing that it strongly overlaps row 12. 

The neighbor sets of a in the figure are {1}, {2}, {3}, {4, 5, ... , 10}. Because it is a P node 
(degenerate), every union of at least one and fewer than all of these sets is a member of M{TZ), 
and this is easily verified using similar checks. 

Only unions of neighbor sets of a single internal node can be members ofJ\f{TZ). For example, 
{1, 2, 10} is not a union of neighbor sets of a single node, and it is not a member of J\f{lZ) because 
it strongly overlaps rows 2, 4, 8, 10, 11, and 14. 

Note that IZ has the circular-ones property, and the figure depicts a way to cyclically order 
the edges incident to each node such that in the resulting tree, every member of TZ is consecutive 
in the circular ordering of leaves. This is an example of a general phenomenon: 

Theorem 2.3 ([IZ]). Let T he the bipartitive tree for J\f{lZ), where TZ is the set family repre- 
sented by rows of a circular-ones matrix with column set C. Then: 

• The edges incident to internal nodes can be cyclically ordered in such a way that the 
resulting cyclic order of leaves is a circular-ones ordering of the matrix; 

• Reversing the cyclic order of edges about a prime node imposes a new circular-ones or- 
dering on the leaves. 

• Arbitrarily permuting the cyclic order of edges about a degenerate node imposes a new 
circular-ones ordering on the leaves. 

• All circular-ones orderings of the leaves are obtainable by a sequence of these two opera- 
tions. 

This gives a convenient data structure, the PC tree, for representing all circular-ones order- 
ings of a circular-ones matrix. 

For example, in Figure [H permuting the counterclockwise cyclic order of neighbors about a 
so that it is (1, c, 2, 3) and reversing the cyclic order of neighbors about b so that it is (c, 10, 9, 8, 7) 
imposes a new cyclic leaf order on the tree, (1, 4, 5, 6, 10, 9, 8, 7, 2, 3), which is easily seen to be 
a circular-ones ordering. 

To gain an insight into why this works, consider a circular-ones ordering of columns of a 
circular-ones matrix, let X be a circularly consecutive set of columns, and let R be the set 
represented by some row. If X is removed and reinserted in reverse order, it will disrupt the 
consecutiveness of R if and only if it strongly overlaps R. Since M{1Z) is the family of sets that 
don't strongly overlap any column, they are the sets that can be reversed in the cyclic order 
without disrupting the circular consecutiveness of any row. Each allowed rearrangement of the 
PC tree corresponds to a sequence of such reversals, and each such reversal is allowed by the 
PC tree. 

The PC tree was first described by Shih and Hsu [28]. Its relationship to bipartitive set 
families and the circular-ones orderings of matrices was first described in p!7] . 
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Figure 1: The PC tree for a circular-ones matrix. Nodes b and c are C nodes (double circles) 
and node a is a P node. 



It takes time linear in the size of a circular-ones matrix to build the PC tree even when a 
circular-ones ordering is not given as part of the input [T71[28]. As part of the output, a circular 
ordering of the edges incident to each internal node is given, which imposes a circular-ones 
ordering on the leaves. This gives a representation of all circular-ones orderings of the matrix 
in linear time. 

Henceforth, we will assume for simplicity that every row of a circular-ones matrix has at least 
one 1 and one 0; any row without this property is irrelevant to the circular-ones arrangements. 



3 Quotient labels for the PC tree 

In this section, we give a way to label the nodes of the PC tree of a circular-ones matrix 
with quotients so that the matrix can be reconstructed from the labeling. These results were 
developed in [8], and the scheme is similar to ones developed for PQ trees in pOj [23 1 125]. 

Recall that we assume that every row of the circular-ones matrix M has at least one 1 and 
at least one 0, since a row that is all zeros or all ones has no effect on the circular-ones orderings 
of the matrix. Suppose that M has a circular-ones ordering of columns. Consider the ordered 
tree that results from the PC tree of M by the cyclic ordering about each node that is forced 
by the cyclic leaf order corresponding to the circular-ones ordering. For each row R that has 
at least two I's we perform the following procedure. Root the PC tree at a leaf corresponding 
to a column with a in i?. We maintain the order about each node u in the following way. If 
{wi,W2, ■ ■ ■ ,Wk) is the counterclockwise order about u and Wi becomes the parent of n, then 
(wj+i, Wj+2, ■ ■ ■ jWhTWi, ■ ■ ■ , Wi-i) becomes the linear order of children in the resulting ordered 
rooted tree. Let u be the least common ancestor of the I's in R. Let X be the leaf descendants 
of a child of u. Since X G J\f{TZ), X does not strongly overlap R. Therefore, X is either a 
subset of R or disjoint from R. It follows that R consists of the leaf descendants of two or more 
children of u. Moreover, since the cyclic order of edges about u gives the circular-ones orderings 
of M, R consists of the neighbor sets of u through a consecutive set R' of children of u. Let R' 
be the projection of R on u. Figure [2] illustrates the concept. 

As a special case, if R has only one 1, let c be the column where the 1 occurs, and let u be 
its neighbor. We consider the projection R' of R to be on u, and to consist of n's neighbor set 
{c}. 

When this projection has been performed on all rows, each internal node u has received the 




Figure 2: The projection R' of a row i? of a circular-ones matrix. Each element of the row 
corresponds to a leaf of the PC tree. Rooting the tree at any leaf where the row has a and 
then finding the least common ancestor of the I's gives the node u that the row projects to. 
The projection makes up a row vector in the quotient matrix at u that has one column for 
each neighbor of u, ordered in cyclic order. Because every neighbor set is a member of AA(7^), 
the columns where a row has I's will always be a union of neighbor sets of u. In this case, 
R has I's in {cs, . . . , C13}, and the neighbor sets of u (dashed arcs) are {co,ci,C2}, {03,04}, 
{cs}, {cg, C7, cg}, {cg, cio, . . . , C13}, and {ci4, C15, . . . , 022}- The ones that are subsets of R are 
reachable through neighbors 0^,3,1. The projection R' of i? is a row of u's quotient matrix 
that has I's in columns 0^,3, t, and O's in columns x,y,z. 

projection of zero or more rows of M, and we represent each projection with a row vector in a 
quotient matrix whose columns are neighbors of u. The row has a 1 in column w if the leaves 
reachable through w are I's of R and a if they are O's in R. The result is a matrix whose rows 
are sets of neighbors of u. Note that the quotient at a node may be empty. 

Figure [3] illustrates the quotients for the example of Figure [TJ (The tree has been rooted at 
node c; the motivation for this is explained below.) 

Lemma 3.1. At a P node u with k neighbors, each row of the quotient matrix consists of either 
k — 1 I's and one or k — 1 O's and one 1. 

Proof. Let M be a circular-ones matrix, and let 7^ be the family of sets of columns obtained by 
considering each row to be the bit-vector representation of a set. 

By the definition of a P node, every union of neighbor sets of u is a member of M{TZ). Let 
R' be the projection of some row R to the quotient at u. Suppose R' consists of more than one 
and fewer than k — 1 neighbors of u. Then R consists of more than one and fewer than k — 1 
neighbor sets of u. Let x and y be two neighbors of u, one in R' and one not. Let X and Y 
be the neighbor sets of u reachable through x and y. Then X U y is a union of neighbor sets 
of u that strongly overlaps R. This disqualifies it as a member oi Af{TZ), contradicting the fact 
that every union of at least one and fewer than all neighbor sets of a P node is a member of 
A/'(7^). D 

The quotient at node a of Figure [3] illustrates the phenomenon. The projections of rows 2 
and 4 exclude only neighbor c, the projection of row 5 excludes only neighbor 1, the projection 
of row 3 excludes only neighbor 3, and the projection of row 1 contains only neighbor 1. 

Lemma 3.2. Any cyclic ordering of the quotient at a P node is a circular-ones ordering of it. 
The quotient at a C node has a unique circular-ones ordering, up to rotation and reversal. 
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Figure 3: The quotients at the internal nodes of Figured! The tree has been rooted at its center. 
Rows {1,2,..., 5} project to node a, rows {6, 7, 8, 9} project to node b, and rows {10, 11, ... , 14} 
project to node c. 



Proof. Again, let M be a circular-ones matrix, and let 7^ be the family of sets of columns 
obtained by considering each row to be the bit-vector representation of a set. 

The result follows immediately for P nodes by Lemma 13. 1[ Suppose u is a Q node and let 
k be the number of neighbors of u. Let X be an arbitrary union of more than one and fewer 
than k — 1 neighbor sets of u. By the definition of a Q node, X N{TZ). Therefore, it strongly 
overlaps some row R oi M. Let y be a column not m. X \J R. Root the PC tree at y. \i R 
projects to a proper descendant of u, then R is either a subset of X or disjoint from X, a 
contradiction to the strong overlap. If R projects to a proper ancestor of li, then X C R, also a 
contradiction. If R projects to a node that is neither an ancestor nor a descendant of u, then R 
is disjoint from X, once again a contradiction. Therefore, R projects on u, and the projections 
X' and R' of X and R on u are strongly overlapping sets of neighbors of u. 

We conclude that every set X of more than one and fewer than k — 1 neighbors of u strongly 
overlaps some row R of the quotient at u. Reversing the order of members of X in the cyclic 
ordering of neighbors about u disrupts consecutiveness of R. This implies that the PC tree of 
the quotient matrix at u has a single internal node, a C node. The circular-ones orderings of 
this quotient are unique up to rotation and reversal. D 

As an illustration, it is easily verified in Figure [3] that every ordering of columns of the 
quotients at P nodes a is a circular-ones ordering, and that the quotient at C nodes b and c 
each have only two circular-ones orderings, one of which is depicted and the other of which is 
its cyclic reversal. 

3.1 Computing the quotient-labeled PC tree in time linear in size{M) 

Note that the linear-time PC-tree construction algorithm of ^ITj gives a circular-ones arrange- 
ment of the leaves, but does not label the nodes with their quotient. Let u be the least common 
ancestor of a row R when the tree is rooted at some leaf not in R. Let a node be black if all of 
its leaf descendants are in R. Because all neighbor sets of u are either subsets of R or disjoint 
from it, the black nodes consist of all ancestors of the leaves in R, up to, and possibly including 
u, if all of its children are black. The black children of u are the ones corresponding to the 
projection of R. 

To find the projection of R, we can therefore blacken the leaves that are members of R. 



When we blacken a node we increment a counter in the parent, so that it keeps track of how 
many blackened children it has. When a node's counter is incremented to a value that equals 
its number of children, we blacken it, and increment its parent's counter. It is easily seen by 
induction on the height of a node that a node is blackened by this procedure if and only if it is 
black. 

The procedure halts when no new vertex can be marked. The blackened nodes induce a 
rooted forest, which is a rooted tree if u is black. Every internal node of this forest or tree has 
at least two black children. Therefore, the number of black internal nodes is bounded by the 
number of black leaves, which is the size of R. Since the procedure spends 0(1) time whenever 
it blackens a node, it takes 0(|i?|) time to find the least common ancestor u and the edges to 
neighbor sets that make up i?'s projection. 

The procedure can be simulated without actually rooting the tree at y; it suffices to work 
"inward" from the marked leaves that are members of R, and mark a node when its counter 
is equal to one less than its number of neighbors, since one of its neighbors is implicitly the 
parent. Since the entry in the quotient matrix has a 1 for each member of X, it takes 0(|i?|) 
time to find the quotient representation of each row R. The total time to find all quotients is 
therefore linear. 

3.2 Computing the quotient-labeled PC tree in 0{n) time 

We now show that given a succinct representation of a circular-ones matrix, we may obtain the 
PC tree in 0{n) time. By complementing a row of a binary matrix, we mean changing every 1 
in the row to a and every to a 1. It is obvious that complementing a row of a circular-ones 
matrix yields a circular-ones matrix [13] . The following is given in [T7] for the version of the 
PC tree where the nodes do not have quotient labels. 

Lemma 3.3. The PC tree of a circular-ones matrix is invariant under complementing a row 
of the matrix. 

The lemma is obvious from the observation that N{1Z) does not change when a row is 
complemented . 

Let c be a column of a circular-ones matrix. Complementing all of the rows that have a 1 in 
column c turns c into a column of zeros. In the succinct representation of the matrix, it takes 
0(1) time to complement each such row: if (/, /) represents the columns of the first and last 1 
in the block of I's in the row, then the {l + l, f — 1), modulo the number of columns, represents 
the complement of the row. Since now no row contains c, a circular-ones ordering where c is the 
last row of the matrix is a consecutive- ones ordering on the columns of the matrix, excluding c. 

The PQ tree was developed by Booth and Lueker to represent all consecutive-ones orderings 
of the columns of a consecutive-ones matrix [3]. It is a certain rooting of the PC tree. (The PC 
tree was originally developed to provide an easier way to compute the PQ tree. [28]) 

Lemma 3.4 ([17]). Let M he a circular-ones matrix and c he a column. Let M' he the result of 
complementing all rows of M that have a 1 in column c, and then removing column c. Removing 
leaf c from the PC tree of M and then rooting it at the neighhor of C gives the PQ tree for M' . 

Lemma 3.5. R takes 0{n) time to compute the PC tree of a circular-ones matrix, given a 
succinct representation of it. 

Proof. Given a succinct representation of a circular-ones matrix M, we may remove all rows 
that have a single 1 or a single 0, since these have no effect on the circular-ones orderings of 
the matrix, hence on the PC tree. We may then select the last column c of M , and identify 



the rows with a 1 in column c in 0{n) time. We may complement each of them in 0(1) time, 
as described above, for a total of 0{n) time. We may remove c from the resulting matrix by 
decrementing the record of the number of columns. This takes 0{n) time and gives a succinct 
representation of a consecutive-ones matrix, M'. 

In [26], an 0{n) algorithm is given for finding the PQ tree of a consecutive-ones matrix, 
given a succinct representation. By Lemma [3.4l and the fact that the removed rows, which have 
all O's or all I's, have no effect on the PQ tree of M', we may then add a new leaf corresponding 
to c to the root of the tree produced by this algorithm, and then unroot it to obtain the PC 
tree of M, also in 0{n) time. D 

Note that since the quotient labels of a PC tree have circular-ones orderings of columns, they 
can be expressed in the form of succinct representations also. Doing this for all quotient labels 
causes them to take 0{n) space, since there is one quotient row for each row of the original 
matrix. This raises the question of whether we can find this succinct quotient-labeled PC tree 
in 0{n) time. 

Lemma 3.6. Given a succinct representation of a circular-ones matrix with n rows, it takes 
0{n) time to find the succinct quotient-labeled PC tree. 

Proof. The algorithm proceeds as in the proof of Lemma 13.51 until the PQ tree of the succinct 
representation of the consecutive-ones matrix M' has been computed. We install the succinct 
quotient labels in this PQ tree, as follows. We use Harel and Tarjan's least-common ancestor 
algorithm [T3] to find the least common ancestor of the endpoints of each interval, and the child 
that contains the right endpoint of each interval. Reversing the tree and repeating this gives the 
child containing the left endpoint. This takes 0{n) time in total. These become the beginning 
and ending points for the representation of the row in the quotient at its least common ancestor. 

For every row that was complemented to obtain M' from M, we complement the image of 
the row in the quotient at its least common ancestor, which takes 0(1) time for each of these 
rows. Note that if the image is all but one neighbor tt; of a node u and w is not a leaf, then, 
after the complementation, this projection consists only of w. However, the definition of the 
projection at the beginning of Section 3 dictates that in this case the projection be onto the 
least common ancestor of the row when the tree is rooted at a node that is in the row. Now 
w, not u, is this node, so the the projection must be moved to w, and be changed to consist of 
all neighbors of w other than u. 

For each row that was removed because it has only one 1 or one 0, let d be the leaf of the 
PC tree (column of the matrix) where the 1 or the occurs. The projection is on the neighbor 
w of d, and consists either of d or of all neighbors of w other than d, depending on whether the 
row had a 1 or a in ti. D 

4 Testing isomorphism of circular-ones matrices 

In this section, we give an 0{size{M)) time algorithm, first developed in [8j, to test isomorphism 
of two circular-ones matrices if no circular-ones ordering of the two matrices is given, and an 
0{n) algorithm if a succinct representation is given. 

It is possible to test whether the matrices are both circular-ones matrices in linear time 
and to find a circular-ones ordering if they are [31]. If both of them fail to be circular-ones 
matrices, the input is rejected for failing to meet the precondition. If exactly one of them fails 
to be a circular-ones matrix, they are not isomorphic. Otherwise, from this ordering, we may 
find a succinct representation of the two matrices in linear time. We may then produce the 
quotient-labeled PC trees with a succinct representation of the quotients in 0{n) time. 
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If the numbers of rows that are all O's in the two matrices differ or the number of rows of 
all ones differ, the matrices are not isomorphic. Otherwise, the problem reduces to the question 
of whether the two matrices are isomorphic when these rows are eliminated. Let Mi and M2 
be these two matrices. This allows us to continue under our assumption that every row of the 
matrices has at least one and at least one 1. 

Let us define an isomorphism vr from one quotient-labeled PC tree, T, to another one, T' to 
consist of the following. It must be that whenever u and v are nodes of T, then Tr{u) and tt{v) 
are adjacent in T' if and only if u and v are adjacent in T. This is the standard notion of graph 
isomorphism. It must also satisfy an additional constraint. Each neighbor w of internal node u 
of T corresponds to a column of the quotient matrix at u. Each neighbor w' of node u' of T' is a 
column of the quotient matrix at u'. If vr is an isomorphism of the underlying trees, then let vr^ 
denote the bijection it induces from neighbors of u to neighbors of u' . It must be the case that 
at each internal node u, vr^ is an isomorphism from the quotient matrix at u to the quotient 
matrix at u' . If such a vr exists, we say that the trees are isom,orphic as quotient-labeled PC 
trees; otherwise we say that they are not. This definition precludes mapping a P node to a Q 
node, because the quotient label determines whether a node is a P node or a Q node. 

Lemma 4.1. Two circular-ones matrices are isomorphic if and only if their quotient-labeled 
PC trees are isomorphic. 

Proof. Suppose matrices M and M' are circular-ones orderings of isomorphic circular-ones ma- 
trices. Let TT be an isomorphism from M to M' . After permutation of columns of M by vr, 
M and M' have identical quotient-labeled PC trees, since they are identical multisets of row 
vectors. Since M and M' are both circular-ones orderings of M, tt is one of the permutations 
of columns of M permitted by the quotient-labeled PC tree of M. In other words, the PC tree 
of M' can be obtained from the PC tree of M by the permitted operations. These operations 
define an isomorphism from the quotient-labeled PC tree of M to the quotient-labeled PC tree 
of M'. 

Conversely, suppose M and M' are circular-ones orderings of matrices that have isomorphic 
quotient-labeled PC trees. Use the isomorphism to arrange the PC tree of M so that it is 
identical to that of M'. Since it now represents M' instead of M, it follows that the permutation 
of leaves of the PC tree induced by the isomorphism is an isomorphism from M to M' . D 

Lemma |4. II is the basis of our algorithm: for two circular-ones matrices, we compute their 
quotient-labeled PC trees and test whether they are isomorphic. To test whether they are 
isomorphic, we encode the trees with strings in such a way that they both have the same 
encoding if and only if they are isomorphic. 

4.1 Testing isomorphism of rooted, unordered trees 

Our starting point is an algorithm for testing isomorphism of rooted, unordered trees that is 
given in the textbook [T]. A rooted tree is unordered if the left-to-right order of children of 
its nodes is not specified. Two unordered, rooted trees Ti and T2 are isomorphic if they are 
isomorphic as directed graphs when the edges are oriented from child to parent. 

If the trees do not have the same height, they are not isomorphic. Otherwise, the algorithm 
proceeds by induction by level, from level 0, which is the level of the deepest leaf, to level h, 
which is the level of the root and the height of the trees. At each level, it labels each node u with 
an integer isomorphism- class label e„ such that two nodes at the level have the same label if 
and only if the subtrees rooted at them are isomorphic. At step i, we may assume by induction 
that this has been done for nodes at level i — 1. For all leaves at level i, the isomorphism-class 
label is 0. For the remaining nodes, we may apply the following procedure. Let tu be the tuple 
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of isomorphism-class labels for the children of u, sorted in nondecreasing order. We may then 
sort the non-leaf nodes at level i by lexicographic order of the tuples assigned to them, in order 
to group identical tuples together. If there are k distinct tuples, among the tuples at level i, we 
then assign isomorphism labels from 1 through k to the tuples, and then each node with the 
number of the tuple it was assigned. By induction, this meets the precondition for induction 
step i + 1 . 

Therefore, two trees are isomorphic if their roots receive the same isomorphism class label. 

4.2 Canonical encodings of quotients 

An additional requirement of an isomorphism vr on quotient-labeled PC trees is that when a 
node u of T maps to a node tt{u) in T', then the neighbors of u map to neighbors of tt{u) in a 
way that is an isomorphism of the quotient at u to the quotient at 7r{u). 

This requires that u and tt{u) have isomorphic quotients. In this subsection, we give an 
encoding of the quotient at a node of a quotient-labeled PC tree so that two nodes receive the 
same encoding if and only if their quotients are isomorphic. 

4.2.1 P nodes 

Consider the P node a of Figure [3l An obstacle to an immediate canonical representation is that 
the quotient can be presented in any column order, since all column orders are circular-ones 
orderings. 

By Lemma 13. H every row of the quotient consists either of one of n's neighbors, or of all 
but one of its neighbors. We can encode the quotient by giving, for each neighbor, an ordered 
pair of integers. The first integer is the number of rows that exclude only that neighbor and 
the second integer is number of rows that contain only that neighbor. 

For example, at node a in Figure [3l the tuple generated by neighbor 1 is (1, 1), since the 
projection of one row, row 5, excludes only neighbor 1, and the projection of one row, row 1, 
contains only neighbor 1. The tuple generated at neighbor 2 is (0, 0) since no row contains only 
2 and no row excludes only 2. The tuple generated at neighbor 3 is (1,0), since the projection 
of one row, row 3, excludes only neighbor 3 and no rows contain only neighbor 3. The tuple 
generated at neighbor c is (2, 0), since the projections of rows 2 and 4 exclude only c and no row 
contains only c. Note that the tuple generated at a neighbor is invariant under permutation of 
the cyclic order of neighbors about the P node. 

Listing these tuples in the cyclic order of the nodes at which they are generated yields an 
encoding of the quotient. The next step of the construction is to sort the generated tuples 
lexicographically to obtain a tuple of tuples. In the example of node a of Figure [3l this gives 
((0, 0), (1, 0), (1, 1), (2, 0)). This encoding is unique, since it is the lexicographic minimum of all 
tuples whose elements are the generated tuples. 

Finally, in order to prevent the possibility that a P node and a C node will have the same 
encoding, we prepend a reserved P-node flag, 0, to the list of tuples. In the example if Figure O 
this gives (0,(0,0), (1,0), (1,1), (2,0)). 

This therefore gives a test of isomorphism of quotients at two P nodes by generating their 
canonical representations and testing whether they are equal. 

4.2.2 A canonical encoding of the quotient at a C node 

In generating a canonical encoding of the quotient at a C node n, there are only two cyclic 
orders of the columns that are circular-ones orderings. One is the reverse of the other. This 
leaves us with two obstacles to a canonical representation: which of these two cyclic orders 
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should we choose, and where in the cychc ordering should we begin in developing a tuple to 
represent the quotient? 

We begin by traveling counterclockwise around the cycle, starting at an arbitrary point. At 
each neighbor w, we generate a tuple that lists the lengths of rows whose clockwise-most 1 is 
at w, and list them in nondecr easing order. 

For example, consider the C-node c depicted in Figure [3j At neighbor a, we see that there 
are two projections, of rows 11 and 10, whose clockwise endpoint is at a, and they have lengths 
2 and 3, respectively. Therefore, the tuple generated for neighbor a is (2,3). At neighbor 4, 
there is one projection, of row 12, that has its clockwise endpoint at 4, and it has length 3, so 
the tuple generated for 4 is (3). Similarly, the tuples generated for neighbors 5, 6, and b are (), 
(2,3), and (), respectively. 

Assembling these tuples in clockwise order, we get ((2, 3), (3), (), (2, 3), ()). However, we 
must consider that we made an arbitrary decision in choosing the point on the circle at which 
to start generating the tuples. The effect of different choices is to rotate the resulting tuple of 
tuples. To choose it in a canonical way, we choose the rotation of the generated list of tuples 
that is earliest lexicographically: ((), (2, 3), (), (2, 3), (3)). We also made an arbitrary decision in 
going around the cycle clockwise instead of counterclockwise. Therefore, we repeat the above 
procedure counterclockwise, generating ((2), (3), (3), (2), (3)) if starting at b, and then choose 
the rotation of this that is lexicographically minimum, ((2), (3), (2), (3), (3)). To choose the 
direction of travel in a canonical way, we choose, from these two lists of tuples, the one that is 
earlier in lexicographic order: ((), (2,3), (), (2,3), (3)). 

To avoid any possibility that a C node and a P node can get the same encoding, we prepend 
a reserved C-node flag, 1, to the list, yielding (1, (), (2), (2,3), (), (3)). 

The general algorithm is as follows. First, we order the neighbors of the node u so that 
the quotient has the circular-ones ordering. For each neighbor w in counterclockwise order, 
list the lengths of the rows of the quotient whose clockwise- most 1 is at tf, in nondecreasing 
order. This gives a tuple {li, h, ■ ■ ■ h) for w. The sequence of tuples generated for each neighbor 
w, taken in counterclockwise order, gives a tuple of tuples. Rotate this ordering to get the 
lexicographic minimum through rotation of the tuples. Then repeat the exercise reversing the 
roles of clockwise and counterclockwise to obtain another such set of tuples. From these, select 
the one that is earlier lexicographically. Then prepend the reserved C-node flag 1 to this list. 

If TT is an isomorphism, then the tuples generated at a neighbor w of u is the same as the 
one generated at neighbor Tr{w) of 7r(n). Since the mapping of the cyclic ordering neighbors of 
u to neighbors of n^u) is the cyclic ordering of neighbors of ir{u) or its reverse, the quotient 
at u and ■k{u) are encoded by the same tuple. Conversely, if the quotient at C nodes u and w 
have the same tuple, then since each tuple uniquely encodes cyclic rotations of the quotient, u 
and w have isomorphic quotients. Therefore, isomorphism of two quotients can be tested by 
determining whether they are encoded by the same quotient. 

4.3 Testing isomorphism of quotient-labeled PC trees 

We combine elements of the rooted-tree isomorphism test of Section ST] with the test for isomor- 
phism of quotients of Section 14.21 i^i order to obtain an isomorphism test for quotient-labeled 
PC trees. 

The use of elements of Section [d] requires us to root the PC trees. Conceptually, a rooting 
of a tree may be viewed as an orientation of its edges from child to parent, yielding a directed 
graph. An isomorphism vr from one rooted tree, T, to another, T', is a bijection from nodes 
of T to nodes of T' such that for u,v G ^(T), {7r{u),Tr{v)) is a directed edge in T' if and only 
if {u, v) is a directed edge in T. Once we root two PC trees, we require them to satisfy this 
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condition. We must therefore be careful to root the two trees in an isomorphic way whenever 
they are isomorphic. 

The center of a one- vertex tree is its vertex and the center of a one-edge tree is its edge. 
Otherwise, the center is obtained by deleting its leaves and recursively finding the center of the 
resulting subtree. It consists of a single vertex or a single edge. 

If the center of a PC tree is node, we root it at that node. If it is an edge vw, we subdivide 
the edge with a pseudo-node x and root it at x so that the tree has a node as a root, as we did 
above. In the quotient at v, replace the name of w with x and in the quotient at w, replace the 
name of v with x. Now x can be treated as a P node with an empty quotient. 

Once we have rooted the trees, we proceed by induction on the level i, as in the algorithm 
of Section 14.11 Because we are applying a stronger notion of isomorphism, which observes 
constraints imposed by the quotients, we must redefine what is meant by the isomorphism-class 
label eu assigned to a node u at level i. 

For every node u in T we define a tree Tu, an induced subtree of T, as follows. If u is a leaf 
node in T, let Tu be the one-node tree consisting of u. If « is an internal node that is not the 
root of T, let Tu be the quotient-labeled subtree induced by u, its descendants, and the parent 
w of u. If u is an internal node and the root of T, then T„ = T, and u is the root of Tu- Two 
nodes u and u' at level i are in the same isomorphism class at level i if and only if T„ and T„' 
are isomorphic as quotient-labeled rooted PC trees. 

Because of the inclusion of the parent of u, the neighbors of each internal node of T^ are 
the same in T„ as they are in T. This allows each internal node to retain the same quotient in 
Tu as it has in T. This tree is rooted at the parent of u. 

In order to merge the constraints of Sections 14.11 and 14.21 we prepend the isomorphism class 
label Cu) of a node w at level i — 1 (Section 14. ip to the tuple generated for the node in the 
encoding of the quotient (Section 14. 2|) . 

The preconditions at the beginning of the induction step at level i are the following. Leaves 
at level i are labeled with equivalence-class label 0. If u is a non-leaf at level i, then the parent 
p of It, if it exists, is a node of T„ and labeled with equivalence class label -1. Note that p is a 
leaf of Tu, but no automorphism T to itself will map p to any other leaf of T„, since they are at 
different levels of T„. Therefore, we must give p a different isomorphism class label from other 
leaves of Tu- 

For any other neighbor w oi u, w lies at level i — 1, and, by induction, it is labeled with 
an integer isomorphism class label for level i — \. The isomorphism classes reflect the stronger 
constraints where, if w and w' are two nodes at level i — \ that have the same label, T^ and 
Tuji are isomorphic as rooted quotient-labeled PC trees. If, together in T and T', there are k 
distinct isomorphism equivalence classes for internal nodes at level « — 1, they are labeled with 
integers between 1 and k, where 1 denotes that a vertex is a leaf at level i — 1. For each neighbor 
of u, let Cu, denote the integer label from {—1,0, 1, . . . , A;} assigned to w. 

We now strengthen the inductive step of Section [4. II to make the stronger induction hypoth- 
esis go through for level i. The canonical encoding of the quotient at u given in Section 14.21 
assigns a unique tuple to each neighbor w of n; to this tuple we simply prepend eui to the tuple 
generated for w. The eui labels on neighbors enforce the constraint that u and u' can get the 
same tuple only if there is a bijection vr from neighbors of u to neighbors of u such that T^, and 
T^(u}) &re isomorphic as quotient-labeled trees. The rest of the tuple forces the constraint that 
they can get the same tuple only if there exists such a vr that is also an isomorphism from the 
quotient at u to the quotient at u' , as in Section H^ Conversely, after ordering the tuples in the 
canonical way described in Section \A.2\ it is clearly sufficient for T„ and T„/ to be isomorphic as 
quotient-labeled trees for u and u' to be assigned the same tuple. 

Replacing the tuples with integer codes from 1 to k, where k is is the number of distinct 
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tuples at level i completes the induction step. 

Therefore, T and T' are isomorphic quotient-labeled PC trees if and only if, after rooting 
them at their centers and performing this algorithm, the roots get assigned the same integer 
label. 

4.4 Time bound 

Theorem 4.2. Given the sparse representations of matrices M and M' , it takes 0{size{M)) 
time either to determine that neither is a circular- ones matrix, or else to determine whether 
they are isomorphic. Given succinct representations of two circular-ones matrices, this problem 
takes 0{n) time to solve. 

Proof. It takes 0{size{M)) time to determine whether they have the same number of I's, by 
counting I's in the two matrices in parallel. If so, size{M) = size{M'). If the standard sparse 
representations of the matrices is given, it takes 0{size{M)) time to determine whether they are 
circular-ones matrices. If neither is, the claim is satisfied. If only one is, they are not isomorphic. 
Otherwise, it takes 0{size{M)) time to convert them to the succinct representations. From the 
succinct representations we can compute the two quotient-labeled PC trees, as described above. 
Therefore, by Lemma [4.H it suffices to show that the quotient-labeled PC-tree isomorphism 
algorithm can be implemented to run in 0{n) time. 

Proposition 1: Summing, for every level i, the number of nodes at level i — 1 plus the 
number of rows in quotients at level i gives a number that is 0{n). This just counts the number 
of nodes in the tree plus the number of rows in the quotients. Each row of a matrix projects to 
just one row of a quotient. 

Proposition 2: At level i, the sum of lengths of the tuples is at most proportional to the 
number of nodes at level i — 1 plus the number of rows in quotients at level i. This is because 
one tuple is generated for each neighbor w of u contains an integer equivalence class label e^, 
and an encoding of a set of rows of the quotient at u. The encoding of each row in the quotient 
only appears in one of the tuples for neighbors. 

Proposition 3: At level i, the maximum integer in any tuple is bounded by the number 
of nodes at level i — 1 plus the number of rows in quotients at level i. The integer equivalence 
classes at level i — 1 are assigned consecutive numbers, starting at 2, by sorting the tuples 
lexicographically, and giving the same integer to two consecutive tuples that are identical, and 
giving an integer that is one higher than its predecessor's if it differs from its predecessor. Each 
row of a quotient at level i maps to only one element of a tuple generated at level i. 

Proposition 4: A radix sort of a set of tuples of integers takes time proportional to the sum 
of lengths of the tuples plus the size of the range of integer values occurring in the tuples [B] . 

The tuples at P nodes must be sorted lexicographically. Number the nodes at level i in any 
order to assign them identification numbers, or I.D. numbers. The maximum number label of 
one of these nodes is at most the number of nodes at level i — 1. To each tuple for a child of a 
P node, prepend the I.D. number of the parent. Sort all tuples for P nodes at level i — 1 in a 
single lexicographic sort. Since the I.D. number of the parent is the major sort key, this groups 
all tuples of children of a P node together, in lexicographic order, giving, for each P node, one 
lexicographically sorted list of tuples for its children. By Propositions 2, 3, and 4, the time for 
sorting all lists of tuples for children of P nodes conforms to the measure given in Proposition 
1. 

The order of tuples of children of a Q node are already given by the canonical procedure for 
generating them, as described above. 
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Figure 4: A non-Helly circular-arc graph. 

We must also sort the set of lists of tuples at level i in order to generate the equivalence 
class numbers for the nodes at level i. A list of tuples can be represented by a simple tuple 
of integers by appending a special separator, -2, to each tuple, and then concatenating them. 
This change of representation does not affect the lexicographic order of the lists, but turns them 
from lists of tuples to lists of integers to make it easier to see that they can be radix sorted. 
The addition of the separators increases the range of values by 0(1). By Propositions 2, 3, and 
4, the time for sorting the set of lists of tuples for children of P nodes conforms to the measure 
given in Proposition 1. 

Assigning integer equivalence-class labels to the lexicographically sorted set of lists of tuples 
trivially takes time proportional to the sum of lengths of lists of tuples, which, by Proposition 
2, conforms to the measure of Proposition 1. 

All of these steps conform to the measure of Proposition 1, so, by Proposition 1, they take 
0(n) time over all iterations of the induction step. 

We must also bound the time to choose, from the Ik possible choices of a list of tuples at 
a C node, one that is earliest in lexicographic order. Generate two lists, one for each cyclic 
ordering, starting at an arbitrary node for each. Turn each of the lists from a list of tuples to a 
list of integers, using the separators, as described above. Then apply the linear-time algorithm 
of ^J to find the cyclic rotation of each list that is earliest in lexicographic order. Of these two 
resulting lists, choose the one that is earlier in lexicographic order. D 

5 Helly circular-arc graphs 

Every interval model has the Helly property [lOj. However, unlike interval models, circular-arc 
models may fail to have the Helly property. Figure Ogives a circular-arc model of a graph where 
the arcs that make up a clique, {^4, B^ C}, do not have a common intersection. 

A graph is a Helly circular-arc graph if it admits at least one circular-arc model that has 
the Helly property. It is easily verified that there is no circular-arc model of the graph of 
Figure m where A, i3, and C have a common intersection point. This illustrates that the Helly 
circular- arc graphs are a proper subset of the circular- arc graphs. 

Not every circular-arc model of a Helly circular-arc graph has the Helly property. Removing 
arcs {D, E^ F} from the model of Figure S] yields a model for the complete graph on three 
vertices that does not have the Helly property. However, this graph is a Helly circular-arc 
graph; it is easy to represent this graph with three copies of the same arc, which has the Helly 
property. 

The strategy of our algorithm for finding whether two Helly circular-arc graphs are iso- 
morphic is related to that of Lueker and Booth for finding whether two interval graphs are 
isomorphic [23]; both algorithms use the clique matrices of the graphs. The main new challenge 
involves correctly computing the clique matrix. 
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We consider adjacency- list representations of two graphs G and G' . Recall that we assume 
that the numbers of vertices and edges in both graphs are the same. For each graph, we obtain 
a Helly circular-arc model or determine that none exists, in 0{n + m) time, using an existing 
algorithm for this problem [18]. If exactly one of them is a Helly circular-arc graph, we determine 
that they are non-isomorphic. If both of the input graphs fail to be Helly circular-arc graphs, 
we reject them for failing to meet the precondition, even though they may be isomorphic. 

If both graphs are Helly circular-arc graphs, we use the Helly circular-arc models to find 
succinct representations of circular-ones orderings of their clique matrices. This involves some 
complications not present in the corresponding problem on interval graphs, which we show 
how to get around in 0{n) time, below. Once we have succinct representations of the clique 
matrices, we use the straightforward fact that two graphs are isomorphic if and only if their 
clique matrices are isomorphic (Lemma 15. 2t below). This reduces the problem to isomorphism 
of circular-ones matrices, which we have solved in 0{n) time in Section [4. 31 The total time is 
0{n + m). 

If the two graphs are circular- arc graphs and the inputs are circular-arc models, we use 
the 0{n) algorithm of [18] to find a Helly circular- arc model for each of them or determine 
that it is not a Helly circular-arc graph. We then proceed as in the case where adjacency-list 
representations of two graphs are given, but take a total of 0{n) time, rather than 0{n + m). 

The main result of this section is the following theorem: 

Theorem 5.1. Given the adjacency-list representations of two graphs, G and G' , where G has 
n vertices and m edges, it takes 0{n + m) time either to determine that neither is a Helly 
circular-arc graph, or else to determine whether they are isomorphic. Given circular-arc models 
of two circular-arc graphs, this problem takes 0{n) time to solve. 

If n' and m' are the number of vertices and edges of G", it takes 0{n-\-m) time to determine 
whether n = n' and m = m', by counting these elements in the two graphs in parallel. If this is 
not the case, then they are not isomorphic. Therefore, we may assume henceforth that n = n' 
and m = m' . 

Lemma 5.2 {^). Two graphs are isomorphic if and only if their clique matrices are isomorphic. 

Proof. A graph isomorphism maps maximal cliques to maximal cliques, so it defines an isomor- 
phism of their clique matrices. Conversely, since two vertices are adjacent if and only if they 
are in a common clique, an isomorphism from the clique matrix of one graph to that of another 
defines a graph isomorphism. D 

In a circular-arc model of a graph, let an intersection segment be a place where a counter- 
clockwise endpoint of an arc is followed immediately by a clockwise endpoint of an arc in the 
model in the clockwise direction; the intersection segment is the region between the two points. 
In an interval model, a set of arcs corresponds to a maximal clique if and only if it is the set 
of intervals containing an intersection segment. Each intersection segment can be located and 
marked by a clique point lying in the segment. A clique is a maximal clique if and only if it 
is the set of arcs that contain a clique point. A consecutive-ones ordering of the clique matrix 
can be obtained by making one column for each such maximal clique, and putting a 1 in the 
column in each row corresponding to an arc that passes through the region. 

By analogy, in a Helly circular-arc model, we can place a clique point in each intersection 
segment. Because the model has the Helly property, a maximal clique must be the set of arcs 
containing one of the clique points. However, in contrast to the case of interval models, not 
all such sets of arcs must be maximal cliques. Moreover, the same clique may appear multiple 
times, as the arcs containing two different clique points in widely separated parts of the circle. 
Figure [5] gives an example. 
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Figure 5: A Helly model where not every intersection segment corresponds to a maximal clique. 
The sets of arcs that contain the intersection segments at points D and F are subsets of the one 
that contain the intersection segment at A. 

To get a clique matrix for the graph, we must eliminate redundant clique points and clique 
points that do not represent maximal cliques from this Helly circular-arc model. Suppose we 
accomplish this. We can obtain a succinct representation of a circular-ones ordering of its clique 
matrix as follows. Number the arcs and label each arc's two endpoints with its arc number. 
Create an array of n buckets, one for each vertex. Number the cliques by numbering the clique 
points that have not been eliminated in order around the circle. This is a circular-ones ordering 
of the clique matrix. 

It suffices, for each bucket i, to store the counterclockwise-most and clockwise-most clique 
number of each arc in order to get the succinct representation of the circular-ones ordering of 
the clique matrix. The counterclockwise points that correspond to clique j are those in the 
maximal consecutive block of counterclockwise points that lie immediately counterclockwise 
from clique point j. The clockwise points that correspond to clique j are those in the maximal 
consecutive block of clockwise endpoints immediately clockwise from clique point j. This can 
be recorded in each bucket, giving the succinct representation of a circular-ones ordering of the 
clique matrix. 

We now describe how to eliminate redundant clique points and clique points that do not 
correspond to maximal cliques. Since we have placed one clique point at each intersection 
segment, and each intersection segment contains the counterclockwise point of an arc and a 
clockwise point of an arc, we have placed at most n clique points. A preliminary version of this 
procedure appeared in [21]. 

For any point p on the circle of a circular-arc model, denote by C{p) the family of arcs of 
the model that contain p. Given two points pi and p2 on the circle of a circular-arc model, let 
us say pi dominates p2 if C{p2) C C{pi). Let A he a circular-arc model and P = {pi, . . . ,pk} 
be a set of points of the circle on which A resides, where (pi, . . . ,pk) is the clockwise order in 
which P appears in a traversal of the circle, starting at an arbitrary point on the circle. Let us 
say that P' C P is a P-dominating set if every point m P \P' is dominated by some point in 
P' . Any minimal set of dominating points, with respect to containment, among the set P of at 
most n clique points we have placed on the circle, is a non-redundant set of clique points. We 
solve the following more general problem in 0{n + |P|) time: 

• Given a set P of points on the circle of a (not necessarily Helly) circular-arc model, find 
a minimal P-dominating set. 

If the model is Helly and P consists of one point per intersection segment, a minimal P- 
dominating set gives the columns of the clique matrix. 

The ascending semi- dominating sequence of P is the subset SD^{P) = {pi E P | C{pi) % 
C{pj) for all pj G P such that 1 < i < j < k}. In other words, SD^{P) contains the points 
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Pi £ P that are not dominated by any later point in P. Similarly, the descending semi- 
dominating sequence of P is the subset SD~{P) = {pj € P \ C{pj) ^ C{pi) for all pi & P 
such that 1 < i < j < k}. The following lemma reduces the problem of finding a minimal 
P-dominating sequence to that of finding SD^ and SD~ . 

Lemma 5.3. Let A be a circular-arc model and P be a set of points on it. Both SD^{SD^{P)) 
and SD^{SD~{P)) are minimal P-dominating sequences. 

Proof. We only prove that P* = SD^{SD^{P)) is a minimal P-dominating sequence. The 
proof for SD~{SD~^{P)) can be obtained by taking the reverse of A. Let P = {pi, . . . ,pk} be 
points on the circle where (pi, . . . ,pk) is the order in which P appears in a clockwise traversal 
of the circle. We first prove that P* is in fact a P-dominating sequence. 

By definition, every point pj G P \ SD~{P) is dominated by some point pi G P for some 
^ ^ i < j- If i is the minimum element in {1, 2, . . . , fe} such that pi dominates pj, then no 
point p € {pi, . . . ,Pi-i} can dominate pi] otherwise p would dominate pj, contradicting the 
minimality of i. Therefore, every point in P\SD~{P) is dominated by some point in SD~{P). 
We can apply a symmetric arguments for SD~{P) and P* to conclude that every point in 
SD~{P)\P* is dominated by some point in P*. Since domination is a transitive relation, every 
point of P is also dominated by some point in P*, i.e., P* is a P-dominating sequence. 

We now show that P* is minimal. Observe that it is enough to show that if a point pi is 
dominated by a point pj € P* , where pj ^ pi, then pi ^ P*. This will imply that no point of P* 
dominates any other point of P*, so no point of P* can be removed from it to yield a smaller 
dominating set. li j < i then pi SD^{P), hence pi P*. If j > i, then since pj G P*, it 
follows that Pj G SD-{P), and we obtain again that pi SD+{SD-{P)) = P*. D 

The algorithms for finding SD^ and for finding SD^ are symmetric. We describe the one to 
find SD^ . The algorithm works by induction on i to find SD~^(Pi), where Pi = {pi,P2, ■ ■ ■ ,Pi}- 
That is, we find those points of Pj that are not dominated by any later points of Pj. 

By induction, assume we have the following partition of SD~^(Pi) at the end of step i: 

• Di'. points in SD~^{Pi) that are already known to be in SD^(P). 

• Qf. remaining points in SD^{Pi); these are points that are not dominated by any later 
point in Pj, but that may or may not be dominated by points in {pj+i, pj+2, . . . ,Pk}- 

It is easy to see that it will follow that when i = k,we get that SD^{Pk) = SD{P) = DkUQk, 
and this solves the problem. 

We begin with SD^{Pi) = {pi}, where Di = and Qi = {pi}. 

In step i + 1, we obtain Qj+i and -Dj+i from Qi and Di as follows. We remove points from 
Qi and insert them in Di if they pass a test that shows that they cannot be dominated by any 
later point in P, including Pi+i. The addition of these points to Di gives -Dj+i. We discard 
other points from Qi that are dominated by Pi+i- We then add pi+i to Qi. This gives Qj+i. 

The test of whether a point q moves from Qi to -Dj+i consists of determining whether it is 
contained in the arc Pj+i that does not contain p^, has its clockwise endpoint in \pi,pi+i), and 
among all such arcs, extends farthest in the counterclockwise direction. (See Figure (H) If q 
is contained in -Bj+i, then, since -Bj+i does not contain any point from {pj_|_i,pj_|_2, . . . ,Pk}, Q 
cannot be dominated by any point in this set, and since q E SD^(Pi), it it is not dominated by 
any later point in Pj. Therefore, it is a member of SD'^{P), and can be moved to -Dj+i. 

We implement Qi as a stack (qi,q2, ■ ■ ■ ,qj), where qj is the top of the stack. Note that the 
set of points that get moved to Di are consecutive at the top of the stack. To move them, we 
pop the stack until we reach an element not in Pj+i, and move the popped elements to Z)j. 
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Figure 6: Computing Dj+i and Qj+i from Di and Qj. Out of all arcs that have their clockwise 
endpoints in [pj,pj_(_i) and do not contain pk, Bi is the one that extends farthest counter- 
clockwise. Elements of Qi that are contained in -Bj+i cannot be dominated by any point in 
{pi+i,Pj4.2, . . . ,Pfc}) hence they are moved to Dj, yielding Dj+i. Out of all arcs that contain p^ 
but do not contain Pi+i, Aj+i is the one that extends farthest clockwise. Elements still in Qi 
that are not contained in Aij^i are dominated by Pi+i, hence discarded from Qi. Then Pi+i is 
added to what remains of Qi, yielding Qj+i. 

A point q' that is still in Qi fails to be dominated by pi+i if and only if it is contained in 
some arc that does not contain pi+i- All arcs in this set contain p^, since otherwise, q' would 
already be identified as a member of Dj for some j < i + 1. 

Of all arcs that exclude Pi+i but contain pk, let Aj_|_i be the one whose clockwise endpoint 
extends farthest clockwise from pk- (See Figure[6j) Since Aj+i is the arc in this set that covers 
the most members of Qi \ -Dj+i, it follows that the points of Qi that are dominated by Pi+i are 
those that are not contained in Aj+i. 

Note that the ones that are dominated by Pi+i are again consecutive at the top of the stack, 
so we pop the stack until we reach an element that is contained in Aj+i, and discard the popped 
elements. 

Since Pi+i belongs in Qi+i, we obtain our stack for Qj+i by pushing pj+i to what remains 
of the stack for Qi . 

For the time bound, note that we may find -Bj+i for each i G {1,2,... A;} by traversing 
[Pj,Pj_i_i) comparing the counterclockwise endpoints of arcs whose clockwise endpoints are in 
this interval and do not contain p]f_. Over all i, this expends 0(1) time on each arc, so it takes 
0{ri) time. 

To find Aj+i for all i € {1,2, ...fc}, we start just counterclockwise from pk and traverse 
the circle counterclockwise, keeping track of the best arc so far. The best arc is initially null. 
When we reach an arc A, we check whether A contains pk, and, if so, whether it extends farther 
clockwise than the best arc so far. If so, A becomes the best arc so far. Each time we reach a 
point Pi, we record the best arc so far as ^i+i. Over all i, this also expends 0(1) time on each 
arc, for a total of 0{n) time. 

The management of the stack implementing Qi takes 0{n) time over all steps, since each 
point is pushed once, and when points are popped, they are consecutive at the top of the stack. 

6 Isomorphism of F circular-arc graphs, convex-round graphs, 
and proper circular-arc graphs 

In this section we show that the circular-ones matrix isomorphism test of Section |4] can be used 
to test isomorphism of F circular-arc graphs and convex-round graphs, using results of Chen [5]. 
From this we get a new algorithm for isomorphism of proper circular-arc graphs. 
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Theorem 6.1. Given adjacency lists of two graphs, it takes 0{n + m) time to either determine 
that the graph are not T circular-arc graphs or to determine whether they are isomorphic. 

Theorem 6.2. Given adjacency lists of two graphs, it takes 0{n + m) time to either determine 
that the graph are not convex-round graphs or to determine whether they are isomorphic. 

A graph is a F circular- arc graphs if its augmented adjacency matrix has the circular-ones 
property. Chen [5J showed that two T circular-arc graphs are isomorphic if and only if their 
augmented adjacency matrices are isomorphic. 

Given adjacency- list representations of two graphs, it takes 0{n + m) time to determine 
whether their augmented adjacency matrices have the circular-ones property \\.2>\ WJ\ [3lJ. If 
they both do, then it takes 0{n-\-m) time, using the isomorphism test of Section H] to determine 
whether they are isomorphic. If a circular-one ordering of the adjacency matrices of the two 
graphs are given using succinct representations, the test takes 0{n) time. 

Convex-round graphs are complements of F circular-arc graphs [2]. The adjacency matrix of 
a convex-round graph has the circular-ones property. Chen [5] showed that two convex-round 
graph are isomorphic if and only if their adjacency matrices are isomorphic. Therefore, we use 
the same technique to get the same bounds for testing isomorphism of convex-round graphs as 
we do for testing isomorphism of F circular-arc graphs. 

Since every proper circular-arc graph is a F circular-arc graph fST], we can use the same 
algorithm also for an isomorphism test of proper circular-arc graphs. This gives a new 0{n-\-m) 
isomorphism algorithm for proper circular-arc graph. 

The 0{n + m) time bound is optimal if the input graphs are given by adjacency lists, but it 
is not optimal if they are given by proper circular-arc models. The algorithm of Lin et al. [22 
for the problem solves it in 0{n) time if the circular-arc models are given. We show how to 
achieve an algorithm with the same time bound. We need to find succinct representations of 
circular-one arrangements of the augmented adjacency matrix of the given proper circular-arc 
models. 

Let Ahe a circular-arc model of G. If A is not a proper circular-arc model, then we can 
convert it in 0{n) time to such a model [27j. The model A can be changed in 0{n) time, 
such that no two arcs cover the circle together [19^ I24j . and the model remains proper. After 
changing the model this way, we index the vertices of G according to the clockwise order 
of their counterclockwise endpoints, starting at an arbitrary endpoint. Since there are no 
arc containment or pair of arcs that cover the circle in A, this indexing gives a circular-ones 
arrangement of the augmented adjacency matrix of G. 

To find the last 1 entry in the row of each vertex of G according to this indexing, in 0{n) 
time, we go clockwise around the circle once, starting at an arbitrary counterclockwise endpoint, 
which belongs to a vertex v. If the next endpoint we encounter is a clockwise endpoint of a 
vertex n, then the last 1 in the row of u is in the column of v. If the next endpoint we encounter 
is a counterclockwise endpoint of a vertex u, we set v = u. We end this traversal when we get 
back to the start point. The first 1 of every row of a vertex in G is found symmetrically. 

We conclude that it takes 0{n) time to find a succinct representation of a circular-ones 
ordering of the augmented adjacency matrix of a graph from a proper circular-arc model of the 
same graph. Once we have succinct representation of the augmented adjacency matrices of G 
and G' we can test them for isomorphism in 0{n) time. 

Theorem 6.3. Given adjacency lists of two graphs, it takes 0{n-\-m) time to either determine 
that the graph are not proper circular-arc graphs or to determine whether they are isomorphic. 
Given two circular-arc models of circular-arc graphs, the same task takes 0{n) time. 
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Figure 7: Two circular- arc graphs G and G' . It is easy to see that the two graphs are not 
isomorphic to each other, since they have different number of edges. 





Figure 8: Normahzed circular-arc models of G and G' . The two models share seven arcs in 
common, and the eighth arc [hold) is flipped between the two models. 



7 Hsu's algorithm for circular-arc graphs isomorphism 

In this section we give an example of two circular-arc graphs that are not isomorphic, but the 
algorithm of Hsu p!H] determines that they are. We begin with few definitions from this paper. 
Let G be a circular-arc graph, without universal vertices, and without any pair of vertices v 
and u such that N[v\ = N[v\. A normalized model of G is a circular-arc model of the graph, 
such that for every two arcs v and u: (1) if N[u\ C N[v\ then the arc of v contains the arc of n; 
(2) if every w € y \ A^fv] satisfies N[w] C N[u\ and every w' € V \N[u\ satisfies A^[?i;'] C N[v\ 
then the arcs of u and v cover the circle together. Every circular graph without universal arcs 
and without any pair of vertices with the same neighborhood has a normalized model [15] . 

Let ^ be a normalized model of G. We get the associated chord model of A by replacing 
the arcs of A with chords. Let A' be an associated chord model of G. The chord model A' 
represents a circle graph Gc, whose vertex set is the same vertex set as of G, and two vertices 
are adjacent if and only if their chords in A' intersect. Although there might be more than one 
unique associated chord model for G, the graph Gc is unique. Hsu [15] defined a type of chord 
model called conformal model for Gc- The chord model A! is a conformal model of Gc- Note 
that we do not repeat the definition of conformal model here, we just give an example for one 
such model. We do not require the definition for our purposes. 

The origin of the mistake in Hsu's algorithm is the statement ''To test the isomorphism 
between two such circular-arc graphs G and G' , it suffices to test whether there exist isomorphic 
conformal models for Gc and G'(j " [l_5i Section 9] , where "such circular-arc graphs" refer to 
circular-arc graphs for which both Gc and its complement are connected. 

Consider the graphs G and G' in Figure [71 It is easy to see that G and G' are not isomorphic, 
since they have different number of edges. Normalized circular-arc models of both graphs are 
given in Figure [8l From the circular-arc models, we can see that the chord model in Figure [9] 
is an associated chord model for both graphs, and hence a conformal model of both Gc and 
G'f^- The graphs Gc = G'q are connected and so are their complements. We conclude that the 
statement above is wrong, and the algorithm of [151 falsely finds that G and G' are isomorphic. 

Hsu J16j noted that the isomorphism of conformal models of Gc and G^ , which the algorithm 
of [15] produces, does give a mapping between vertices of G and G' , if Gc is inseparable 
with respect to modular decomposition. However, we do not know how to handle the case 
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Figure 9: The associated chord model of both G and G' . 

when this condition is not satisfied. We note that the algorithm of [15j works correctly for 
isomorphism of circle graphs. With the 0{n'^) time recognition algorithm for circle graphs |30j . 
the isomorphism test takes 0{n'^) time if the graphs are given as adjacency matrices. The recent 
0{{n + m)a{n -\- m)) circle- graph recognition algorithm [T2] leads to the same running time for 
circle-graph isomorphism, where a(-) is the inverse Ackermann function. If chord models are 
given as an input, then the running time of the isomorphism test can be reduced to 0{n + m) 
using techniques similar to those used in p2j and in our paper. 
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